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We study the interplay of quantum statistics, strong interactions and finite temperatures in the two-component 
(spinor) Bose gas with repulsive delta-function interactions in one dimension. Using the Thermodynamic Bethe 
Ansatz, we obtain the equation of state, population densities and local density correlation numerically as a func- 
tion of all physical parameters (interaction, temperature and chemical potentials), quantifying the full crossover 
between low-temperature ferromagnetic and high-temperature unpolarized regimes. In contrast to the single- 
component, Lieb-Liniger gas, nonmonotonic behaviour of the local density correlation as a function of temper- 
ature is observed. 



The experimental realization of interacting quantum sys- 
tems using cold atoms has reignited interest in many-body 
physics of strongly-interacting quantum systems in and out 
of equilibrium Effectively one-dimensional realizations 
of bosonic ^^Rb quantum gases with tunable local interac- 
tion strength 121 [3] IH E] |6l realize the single-component 
Lieb-Liniger model fT', ISl, for which the crossover from 
weakly- to strongly-interacting physics is experimentally ac- 
cessible and well understood from first principles. Observed 
finite temperature thermodynamics 19J even fit predictions 
from the Thermodynamic Bethe Ansatz (TBA) ifTOl . For 
a single bosonic species in one dimension, statistics and 
interactions are intimately related: the limit of infinitely 
strong interactions (impenetrable bosons) causes a crossover 
from bosonic to effectively fermionic behaviour fTP 12^ for 
density-dependent quantities. Density profiles and fluctua- 
tions accessible from exact thermodynamics allow to discrim- 
inate between these fermionized and quasicondensate regimes 
ifTSi fT4l [TSl [161 . Multicomponent (spinor) systems however 
provide a much larger number of different regimes than their 
single-component counterparts, and realize situations where 
important interaction and quantum statistics effects coexist 
and compete. Their thermodynamics has not been extensively 
studied using exact methods; in this work, we wish to high- 
light some unexpected features inherent to a system in this 
class. 

The experimentally realizable ifTTll case of two-component 
bosons in ID with symmetric interactions, which we will fo- 
cus on here, contrasts with the Lieb-Liniger case in many 
ways. The ground state is (pseudo-spin) polarized flS |T9| 
(ferromagnetic), as expected from a general theorem valid 
when spin-dependent forces are absent |(20|, and thus coin- 
cides with the Lieb-Liniger ground state. On the other hand, 
excitations carry many additional branches, starting from the 
simplest spin-wave-like one. These excitations are difficult 
to describe in general, even in the strongly-interacting limit 
(there, no effective fermionization can be used, since the 
two pseudospin components remain strongly coupled), where 
spin-charge separation occurs ETl l22l l23l l24l l25l . The ther- 
modynamic properties are drastically different from those of 
the one-component Lieb-Liniger gas |26| and at large cou- 
pling and low temperature correspond to those in an isotropic 
XXX ferromagnetic chain |27J. Temperature suppresses the 



entropically disfavoured polarized state, and opens up the 
possibility of balancing entropy and quantum statistics gains 
(from the wavefunction symmetrization requirements) with 
interaction and kinetic energy costs in the free energy. Us- 
ing a method based on the integrability of the system, we find 
that this thermally-driven interplay leads to a correlated state 
with many interesting features, the most remarkable being a 
nonmonotonic dependence of the local density fluctuations of 
the system with respect to temperature or relative chemical 
potential. 

For definiteness, we consider a system of N particles on a 
ring of length L, subjected to the Hamiltonian 



Hn — 



N 



92 



2m ^ dx^ 



9lD ^ H'Xt 
l<i<j^N 



Xj). (1) 



The effective one-dimensional coupling parameter giu is re- 
lated to the effective ID scattering length aiD 1281 via the 
relation giD — h^aiD/2m, and to the effective interaction 
parameter -f — c/n (where n = N/L is the total linear 
density) via c — gium/h^. We set h = 2m = 1 to sim- 
plify the notations. Yang and Sutherland ll29l l30ll showed 
that the repulsive delta-function interaction problem admits 
an exact solution iiTespective of the symmetry imposed on the 
wavefunction, meaning that the wavefunctions of ([T]i are of 
Bethe Ansatz form whether the particles are distinguishable, 
or mixtures of various bosonic and fermionic species |31|. 
The ground state and excitations of multicomponent fermionic 
system were studied, both for repulsive and attractive interac- 
tions, by Schlottmann 1321 [33l . but these results cannot be 
translated to the bosonic case we are interested in. 

Specifically, specializing to N atoms of which M have (in 
the adopted cataloguing) spin down, the Bethe Ansatz pro- 
vides eigenfunctions fully characterized by sets of rapidities 
(quasi-momenta) kj, j = 1,...,N and pseudospin rapidities 
Xa, a — 1, M, provided these obey the N + M coupled 
Bethe equations 
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for j = 1,...,N and a = 1,...,M. For a generic eigen- 
state, the solution to the Bethe equations is rather involved. In 
general, the kj rapidities live on the real axis; the are on 
the other hand generically complex, but arranged into regular 
string patterns, each type of string representing a distinct type 
of quasiparticle. The spectrum of the theory thus contains in- 
finitely many branches, with an increasing effective mass. In 
the continuum limit N oo, L oo, N/L fixed, the TBA 
allows to exploit the condition of thermal equilibrium to ob- 
tain the Gibbs free energy of the system as a function of the 
temperature T and of the total /i (= ^^"^^^ with /i^ the chem- 
ical potential specific to the zth component) and relative 
{= ) chemical potentials. As detailed for example in 

ll34l the Gibbs free energy density is given by 
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where e(A) is dressed energy, which is self-consistently cou- 
pled to the length-n string dressed energy e„(A), n = 1,2, ... 
via the system 
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with the standard convolution notation g * h{X) = 
J^^dX'giX - A')/i(A'), and the kernels a„(A) = 

¥ (»cff+A^ and /(A) = ,J,l^l/,y This set of coupled 
equations is complemented with the asymptotic condition 

lim„^oo ^^4^ ~ ^'^ view of its nonlinear nature, we 
solve the infinite system of coupled integral equations numer- 
ically. In order to validate our results, we have in fact imple- 
mented two independent algorithms, one based on fast Fourier 
transforms and the other based on adaptive-lattice integration 
specifically tailored to the computation of the free energy or 
of its various derivatives. While extremely challenging, it re- 
mains possible to obtain good accuracy throughout parameter 
space, except in the extreme limit of vanishing relative chem- 
ical potential for il << T, /i, c or the low nonzero tempera- 
ture limit < T << il,/i,c|35|. Leaving the details to a 
future publication, we here simply concentrate on the results 
obtained. 

We focus on the polarization, defined as (711 — 712) /{rii + 
712) where Ui = Ni/L the linear density of the ith bo- 
son component, and on the local density-density correlator 

' — — ,^ — . The linear densities themselves, to- 

gether with other equilibrium quantities such as the entropy, 
will be discussed more extensively in a future publication. 
Fig. [Tjcontrasts two single-component gases living in two dif- 
ferent traps (Lieb-Liniger gases) against cohabitation in the 
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FIG. 1: Population densities (top) and polarization (bottom) of the 
spinor Bose gas as a function of temperature, for fixed chemical 
potentials, and contrasted to separate Lieb-Liniger gases at corre- 
sponding chemical potentials. The Lieb-Liniger result for the major- 
ity chemical potential is recovered only at zero temperature, when 
ferromagnetism causes complete polarization. 
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FIG. 2: Polarization as a function of the relative chemical poten- 
tial n, (top) at fixed temperature and for different values of the total 
chemical potential fj, and (bottom) at fixed total chemical potential fi 
and for different temperatures. 



same trap, for a generic choice of chemical potentials. The 
polarization of the ground state is clearly visible and gets sup- 
pressed with temperature at a rate depending on /i, which it- 
self sets the effective coupling 7. In the limit of /i ^ 0, 7 
becomes bigger than l/fc^T, and the gas becomes paramag- 
netic. For ^ ^ ( 7 <C 1), the gas shows ferromagnetic 
behavior. The effect of the relative chemical potential on the 
polarization is shown in Fig. [2] Due to the finite temperature, 
the limit il ^ is always unpolarized. Fig. [3] shows data 
for a wide range of the total chemical potential. For /i small 
enough (7^1, paramagnetic), the value of the polarization 
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FIG. 3: Polarization as a function of the total chemical potential ^, 
for fixed relative chemical potential Q. and for different temperatures 
(main) and (inset) for fixed temperature and for different values of 
the relative chemical potential f2. Lowering the temperature or in- 
creasing SI increases the polarization at any ^, leading back to Lieb- 
Liniger behaviour. 



depends only on Q. and T. In contrast, for /i ^ 0, the system 
behaves as a quasi-condensate (7 ^ 1, ferromagnetic). 

The observable g^'^\ which can be obtained from the inter- 
action parameter derivative of ([3]), quantifies the fluctuations 
of density at a local point of space and time. Fig. |4]presents 
data for this quantity as a function of the effective interaction. 
In the limit 7 — > and in the decoherent regime (where the 
reduced temperature r = (^^^'^^^■^2 ^ 1), the value saturates 

between 2 (for 51 00) and 1 + (for ^ 0, where 
A^c = 2 is the number of components), generalizing the Lieb- 
Liniger result |15]. Our data fit well with this prediction, as 
seen in Fig. |4j where for 7 ^ 10~^, g^'^^ approaches 1.5 for 
il and /i small. For bigger /i, the data does not reach this 
value since the reduced temperature is too low and decoherent 
regime is thus not yet reached. On the other hand, when the 
gas becomes strongly interacting, g^^^ vanishes as expected. 

When studied as a function of temperature, the behaviour 
of g^^) is markedly different in the two-component gas than 
in the single-component case. For fixed chemical potentials, 
it can remarkably exhibit a maximum at finite temperature as 
shown on Fig. |5] In the regime of large relative chemical po- 
tential, the gas is polarized and the correlation is monotonic in 
temperature; however, one can clearly observe that for — > 
a peak appears. For this range of temperatures, r is > 1 but 7 
passes from high values at low T to almost zero for high tem- 
perature. As a function of the relative chemical potential (Fig. 
|6| ), exibits a maximum followed by a local minimum. 
The gas is in the decoherent regime at low relative chemical 
potential and the saturation value of the correlation for this 
regime increases with Q.. For bigger values 1 ^ 7 ^ r 
and for > fi, the first component is quasi-condensating and 
g(^) ^ 1. 

The existence of a maximum and minimum of the density 
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FIG. 4: Local pair correlation g^^'' of the spinor Bose gas as a func- 
tion of the effective coupling 7, at fixed temperature and for three 
different values of the total chemical potential /j.. The relative chem- 
ical potential Q is set to a low value. Inset: the same curves for 
the Lieb-Liniger gas. The asymptotes 7 — > differ, but the general 
shape of the curve is very similar. 
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FIG. 5: Top: the local pair correlation g^^^ as a function of tem- 
perature, for fixed total chemical potential g, and four different val- 
ues of the relative chemical potential Q. The non-monotonicity of 
g^^^ in the spinor gas is to be contrasted to its monotonicity in the 
Lieb-Liniger case (top, inset). Bottom: polarization as a function of 
temperature, for the same values of Q and g. At zero temperature, 
polarization is total irrespective of the chemical potential (see also 
Fig. [TJ, illustrating the ferromagnetic-like physics involved in the 
spinor gas. 



fluctuations is the result of an interesting competition between 
interaction energy and entropy. On the one hand, the bosonic 
nature favours ferromagnetic correlations, which in general 
set the small-temperature thermodynamic properties. On the 
other hand, the reduced entropy associated to the polarized 
states and the enhanced spatial density resulting from quasi- 
condensation bear a free energy cost which can or cannot be 
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FIG. 6: Top: the local pair correlation g^^^ as a function of the rel- 
ative chemical potential SI, for fixed total chemical potential fi and 
four different values of the temperature. Bottom: polarization as a 
function of the relative chemical potential, for the same values of T 
and p. 



afforded depending on the value of temperature and of the 
chemical potentials. Our work clarifies and quantifies these 
effects fully throughout the available parameter space: de- 
pending on the precise values of these three parameters, we 
see that the system equilibrates to a state with markedly dif- 
fering correlations. As a corollary, the nonmonotonicity found 
in g'^'^' could point to other interesting consequences for ex- 
perimental realizations of such a system. The population den- 
sities of a two-component Bose gas in a trap, obtainable by 
coupling our method to a local density approximation, would 
also display interesting correlated behaviour as a function of 
position. The ferromagnetic tendencies of the system would 
tend to drive phase separation, leading to an observable corre- 
lated enhancement and depletion of the spatial density profiles 
of the different bosonic species. 
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